Report for Model 6

— for both jointly and separately

2024-08-10

1. Joint Modeling

1.1 Trace plots for convergence check

The current MCMC setting is:

  • 100,000 iteration;
  • 90,000 burn-in;
  • 10 thinning.

1.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum       1.01       1.03
## LDevsum       1.01       1.03
## dh0           1.44       2.21
## dh1           5.72      10.41
## dh2           3.96       7.16
## dl0           1.31       1.75
## dl1           1.37       1.89
## dl2           1.25       1.63
## dl3           1.02       1.04
## 
## Multivariate psrf
## 
## 5.44

1.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

1.4 WAIC results

LevelH LevelL
DIC 1405.5374 24268.361
DIC3 1243.8191 21376.472
PWAIC 89.9428 1524.853
WAIC 1312.6123 22101.590

2. Separate Modeling

2.1 Trace plots for convergence check

The current MCMC setting is:

  • 100,000 iteration;
  • 90,000 burn-in;
  • 10 thinning.

2.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum       1.01       1.02
## LDevsum       1.02       1.06
## dh0           4.37       7.72
## dh1           7.31      12.81
## dh2           3.59       6.71
## dl0           1.42       2.11
## dl1           1.42       2.07
## dl2           1.12       1.33
## dl3           2.24       3.69
## 
## Multivariate psrf
## 
## 7.27

2.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

2.4 WAIC results

LevelH LevelL
DIC 1377.23924 24118.439
DIC3 1243.55301 21367.683
PWAIC 89.58375 1524.755
WAIC 1311.70933 22092.162